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The generic changes of the electronic compressibility in systems which show magnetic instabilities 
is studied. It is shown that, when going into the ordered phase, the compressibility is reduced by an 
amount comparable to the its original value, making charge instabilities also possible. We discuss, 
within this framework, the tendency towards phase separation of the double exchange systems, the 
pyrochlores, and other magnetic materials. 
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I. INTRODUCTION 

The theoretical possibility of phase separation (PS) in 
magnetic systems was first discussed in connexion to the 
Hubbard model for itinerant ferromagnetism and anti- 
ferromagnetism Recent arguments suggest that PS 
is also likely in other magnetic materials, such as those in 
which magnetism is due to double exchange interactions 
Spin polarons, which can be viewed as a manifesta- 
tion of PS on a small scale, have been analyzed in relation 
to the pyrochlores || . Finally, there is an extensive liter- 
ature on PS in two dimensional doped antiferromagnets 
(see, for instance |10-15|), although there is no definitive 
consensus on its existence. 

A variety of different experiments show features which 
are compatible with electronic phase segregation, at least 
on small scales, near magnetic transitions. For instance, 
near an antiferromagnetic transition there is ample evi- 
dence for the formation of charged stripes, both in the 
cuprates [|l6| and in the nickelates jlT]] . In the ferromag- 
netic manganites, various experiments suggest the exis- 
tence of polarons near the Curie temperature p8[ , which, 
as mentioned above, can be viewed as phase separation 
on small scales. In addition, and also in the ferromag- 
netic manganites (i. e. Lai-xCa^MnOs, x w 1/3), inho- 
mogeneous textures |l9| and hysteretic effects [M have 
been reported near T c . Note that we are discussing the 
relatively simple case of the optimally doped, ferromag- 
netic compounds, avoiding the complications which arise 
between the competition between antiferromagnetism, 
charge ordering and double exchange in these materials. 
Finally, there is evidence for inhomogeneous textures in 
ferromagnetic pyrochlores |2jJ |, The purpose of this pa- 
per is to show the theoretical foundations which make 
likely the existence of electronic phase segregation near 
magnetic transitions. We do not pretend to give an ex- 
haustive list of experimental evidence which support this 
view. 

In the following, we analyze how a magnetic transi- 
tion influences the electronic compressibility. A general 



framework is presented in the next section. Simple appli- 
cations to variations of the Hubbard model are presented 
in section III. Then, it is shown that the same framework 
predicts the existence of PS in the pyrochlores (section 
IV). Section V includes in the same framework the known 
results about PS in double exchange models, in the light 
of the present framework. Section VI analyzes the role 
of long range interactions, and some physical quantities 
likely to be affected by a reduced compressibility. Finally, 
the main conclusions of our work are presented in section 
VII. 



II. GENERAL FEATURES OF THE 
ELECTRONIC COMPRESSIBILITY NEAR A 
MAGNETIC PHASE TRANSITION. 

The formation of an ordered phase induces a decrease 
in the free energy of the material, usually called the con- 
densation energy. Because of it, quantities which depend 
on the variation of the free energy with temperature, such 
as the specific heat, show anomalous, non analytic behav- 
ior at the critical temperature. For instance, the specific 
heat shows an abrupt reduction at T c in mean field the- 
ories. The phase transition can also be tuned by varying 
n, the electronic content (per unit cell), in many systems 
leading to phase diagrams like that shown in fig. (|,|,|). 
Thus, one expects an anomalous dependence of the free 
energy on electronic density as the phase boundary is 
crossed by changing the electronic concentration. 

We can get a simple estimate of the effect by using a 
standard Ginzburg-Landau expansion for the free energy: 
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where s is the electronic magnetization and n is the elec- 
tronic density. We now neglect spatial fluctuations, and 
obtain a mean field approximation to T: 
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where we have expanded the dependence of the free en- 
ergy of the paramagnetic phase on n. kq is the (scaled) 
electronic compressibility in the paramagnetic phase, 
and no defines the equilibrium density in the absence 
of magnetism. When T < T c , the magnetization is: 
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and the free energy becomes: 
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Let us now fix the temperature T and expand this ex- 
pression around the density n c such that T — T c (n c ). 
We obtain: 
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And by taking derivatives, we have: 
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The compressibility has a jump at the transition. The 
origin of this discontinuity is the same as that in the 
specific heat. Moreover, it is reasonable to expect that 
this anomaly will be enhanced when fluctuations in the 
critical region are taken into account. 

Alternatively, and keeping T fixed, one can minimize 
eq.(||) first with respect to n. By expanding T c (n) around 
no, we find: 
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and, inserting the value of n in the free energy when 
T<T r : 
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The dependence of T c on density leads to a negative quar- 
tic term in the dependence of the free energy on the 
magnetization. When b/A - (a 2 k )(dT c /dn) 2 /8 < 0, the 
magnetic transition becomes first order. This condition 
is equivalent to saying that the effective compressibility, 
defined in eq.(|J) becomes negative. Thus, PS near the 
transition can be thought of as arising from the transmu- 
tation of a continuous phase transition into a first order 
one by the introduction of an additional field, the density 
n, which is a well known possibility in statistical mechan- 
ics p2| . The new feature found in a magnetic transition 
is that the correction to the compressibility can easily 
be comparable to the initial compressibility. The latter 
is determined by the density of states at the Fermi level 
in the paramagnetic phase. In typical magnetic systems, 
the transition is driven by a coupling constant which is 



of the order of the inverse of the density of states. Fi- 
nally, the dependence of the critical temperature on the 
electronic density depends on the change of the coupling- 
constant with variations in the density of states, which 
is of the same order (see examples below). Thus, no fine 
tuning of parameters is required to obtain corrections to 
the compressibility of the order of the compressibility it- 
self. 



III. PHASE SEPARATION IN HUBBARD LIKE 
MODELS 

A. Ferromagnetic transitions 

We first consider a one band Hubbard model with near- 
est neighbor ferromagnetic exchange couplings: 
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where i is a site index, = (<4 s c i,s)i Si — 

QDs s ' c \ gOs.s'Cj.s') an d <J s ,s' are Pauli spin matrices. 
This model (with appropriate additions) has been in- 
voked as a reasonable starting point for the study of 
itinerant ferromagnetism p3j . Within a mean field ap- 
proximation, the Stoner criterium gives a feromagnetic 
instability for UT>(eF) > 1, where T>(e) is the density of 
states per spin at energy e. Let us assume that in the 
paramagnetic phase there are no electrons per site, with 
chemical potential no- 
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If we shift the two spin bands by ±5, the induced polar- 
ization s — n^i — n^j satisfies: 
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where T>' and T>" stand for the derivatives with respect 
to energy of the density of states. Because of the lack of 
electron- hole symmetry when these derivatives are finite, 
the chemical potential in the polarized state is shifted: 
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The ground state energy, at zero temperature, can be 
written as: 
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where Eq = 2 f e T>(e)de. This is the Ginzburg-Landau 
expansion needed to study the phase transition as func- 
tion of electronic density at zero temperature. The sys- 
tem becomes ferromagnetic when 2?(/i c )[7 > 1. We can 
expand the quadratic term in the magnetization around 
n c = 2 fg c V(e)de as: 
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which leads to the following scaled inverse compressibil- 
ity, as the transition is approached from the ordered side: 
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where k Q 1 = {2V{p, c ))- 1 + (2?7 C ) _1 , and U c = 
(2?(/Lt c )) _1 is the critical coupling for the transi- 
tion to take place. The compresibility is negative 
if V"(fi c )V(n c )/(3V' 2 (n c )) > 1/2. In particular, 
near a saddle point in the 3D dispersion relation, 
we have T>(e) = T>q — c^/]e — egu\. This implies 
lim € ^ esD {V"(e)V(e)/V 2 (e)} -> +oo, and the system is 
always unstable versus PS. 

It is interesting to note that, using a different formal- 
ism, phase separation has been shown to appear near fer- 
romagnetic phases of the Hubbard model in two dimen- 
sions |24|] , in good agreement with the picture presented 
here. 



B. Doped antiferromagnets 

It is well known that the Hubbard model (eq.Q at 
half filling, in a bipartite lattice and with nearest neigh- 
bor hoppings only, has an antiferromagnetic (AF) ground 
state, except in one dimension. The main physical fea- 
tures of this state, a charge gap, long range magnetic or- 
der and low energy spin waves, are well described using 
standard mean field techniques and the Random Phase 
Approximation. Recent work shows that, in the presence 
of a static magnetization, the opening of a charge gap 
occurs while the quasiparticle residues remain finite, at 
least in the infinite dimension limit [ p5| , further support- 
ing the validity of a mean field ansatz. Alternatively, a 
calculation of the ground state energy in the limit of large 
dimensions and small values of U/t can be obtained by 
standard perturbation theory around a mean field sym- 
metry breaking state [^6|. This approach does suggest 
the existence of PS near half filling 10 1. In the follow- 
ing, we analyze the stability of homogeneous mean field 
solutions in arbitrary dimensions, following a different 
approach to that in |hJ . 

The mean field solution at half filling is straightfor- 
ward to extend to finite fillings. Hartree Fock calcula- 
tions show that there are self consistent homegeneously 
doped solutions for a finite range of fil ling s around half 
filling, in two and three dimensions pTipq] . As reported 



earlier, these solutions acquire a negative contribution to 
the electronic compressibility near the transition bound- 
aries, making PS possible. In the following, we analyze 
the stability of homogeneous mean field solutions in arbi- 
trary dimensions, following a different approach to that 

in [nj|. 

The AF distortion shifts the mean field levels, e — > 

sgn(e)-\/e 2 + A 2 , opening a gap and leading to the fol- 
lowing staggered magnetization: 
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where T>(e) = D(— e) and fi are, respectively, the density 
of states (per spin) and Fermi level in the paramagnetic 
phase, with filling factor n — 2 f_ Wg dcD(c), and half- 
bandwidth Wq. The energy is given by 



/ft 
de T>(e) sgn(e) \/ e 2 + A 2 + As 
-Wo 



+ j(n 2 -s 2 ) 
(16) 



with the order parameter satisfying the following self- 
consistency requirement: 
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Both energy and self-consistency equation can be ob- 
tained from the minimization of the following functional 
of the order parameter .F(A), which we take as the start- 
ing point for the stability analysis: 
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Away from half-filling, a critical value U c is re- 
quired for the instability to take place, with U c {n) — 
— 1/ J^ w de -^i. Close to the transition line U c (n), a 
Ginsburg-Landau analysis is straightforward, leading to: 
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The scaled inverse compressibility satisfies: 
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Following the transition line on the AF side to the limit 
of small doping, one has -> 1/T>(fi), leading to 

— > — • Therefore we have proven the existence 
of PS in the vicinity of U c (n), at least close to half filling. 

In fact, PS is not restricted to the previous region. It 
is also present around half-filling, n — > V^, for all values 
of U, as we now show. Close to half-filling and at finite 
U, the compressibility satisfies 
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and, given the following inequality (see eq. 
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proving PS around half-filling for all values of U, as pre- 
viously stated. This analysis indicates that PS is a very 
robust feature of the AF instability at the mean field 
level. In fact, for a flat model density of states, numeri- 
cal results suggests that the entire AF region has negative 
compressibility. 

The study of this section can be extented to the case of 
a diverging density of states at half-filling: 23(e) ~ |e| Q , 
(0 < a < 1). The existence of PS is also proved in this 
case, at least around the transition line U c (n — ► 1). In 
particular, this includes the mean field solution of the 
Hubbard model in the square lattice, where the density 
of states has a van Hove logarithmic divergence. The 
conclusions of our analysis are consistent with the ten- 
dency towards inhomogcncous solutions found in mean 
field studies of the Hubbard model [E9 |32l . 



IV. PYROCHLORES 

The pyrochlores are metallic oxides, MH.2TI2O7, which 
show various anomalous transport and magnetic prop- 
erties, including colossal magnetoresistence. The mag- 
netism is mostly due to the spins of Mn ions. The elec- 
tronic carriers are assumed to come from a wide s band 
from Tl orbitals, with a very low occupation. These 
electrons are coupled ferromagnetically to the Mn spins 
p3|pl| . We can describe the coupled system by the 
hamiltonian 101: 
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where J and J' are positive (ferromagnetic). The model 
undergoes a transition to an ordered where the local spins 
have a finite magnetization, s = | (s) | 7^ 0, and the elec- 
tron gas is polarized, s = |(cj s <7 ss /Cj lS ')| 7^ 0. 

In principle, the entire issue of PS in the py- 
rochlores could be settled using results of Majumdar- 
Littlewood(ML) M. These authors have shown that 
magnetic polarons (i.e. self-trapped, polarized electrons 
in a magnetic bubble of core spins embedded in a param- 
agnetic environment) are the stable configuration of car- 
riers in the paramagnetic phase at low densities. Now 
we will show that the existence of polarons in the para- 
magnetic phase is a sufficient condition for the thermo- 
dynamic instability of the uniform phase of the standard 
mean field analysis. 

Let us consider a low concentration (n) of carriers in 
the polaronic configuration. The free energy at fixed tem- 
perature is proportional to the number of carriers (po- 
larons) : 



F P oi{T,n) = Sf pol (T) 



(28) 



where 8f pol (T) < is the free energy reduction per 
carrier due to its polaronic bubble. The preference for 
polarons implies that the homogeneous (paramagnetic) 
phase Thom is above J- po i in the limit n — > 0. Con- 
sider now that the carrier density reaches the value n at 
which all volume is filled with polarons, each one with 
its confined carrier. It is obvious that, allowing carriers 
to delocalize (keeping core magnetization and tempera- 
ture constant), their kinetic energy diminishes. This im- 
plies the existence of a homogeneous phase (with finite 
core magnetization and total carrier polarization) with 
energy below the polaronic phase at finite carrier con- 
centration. A variational reasoning guarantees that the 
true homogeneous solution of the mean-field treatment is 
even below in free energy. In conclusion, if the single po- 
laron is stable, the free energy of the homogeneous phase 
must satisfy 

F P oi(T, n -> 0) = Sf pol (T) n < T hom {T, n — > 0) 
T P oi{T,h) = Sf pol (T) n > T hom {T,n) 

Remenber that T po i is linear in n, and both T po i and 
Fhom start from the same point. Then, the previous in- 
equalities can only happen if the curvature of T po i with 
concentration changes sign, leading to a negative com- 
pressibility and PS. Therefore we have proven our original 
assertion: stable polarons imply thermodynamic instabil- 
ity towards PS. Notice that this is a sufficient condition, 
but not necessary: PS can exists (see below) even if the 
ML polarons are not the stable configuration in the para- 
magnetic phase. A graphical description of the previous 
argument comparing polaronic and homogeneous phase 
free energies for kT = 0.115 t is shown in fig. [y (two 
upper curves). In addition, we also show the free energy 
of the correct treatment, including carriers entropy (see 
below). 
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FIG. 1. Free energy versus electron concentration for the 
exact mean field treatment (continuous line), the polaronic 
ansatz (dashed straight line), and the degenerate approxima- 
tion for electrons (dashed-dotted line), for kT — 0.125 t, and 
parameters explained in the text. The arrow marks the onset 
of negative curvature. 



There is, however, an important source of concern with 
the previous conclusion. In the ML analysis, carriers in 
the homogeneous phase enter with energy but no en- 
tropy. This would be correct for the usual case of de- 
generate fermions. But the temperature is finite and the 
interesting region corresponds to very low carrier con- 
centration, making the degenerate assumption question- 
able. For the parameters expected to apply in the py- 
rochlores the classical limit would be a more appro- 
priate starting point. In fact, in the limit of zero con- 
centration at finite temperature, the chemical potential 
(n = d n Thom) corresponds to classical particles, and di- 
verges: fi(n — > 0) ~ log(n). This implies that, in that 
limit, the homogeneous paramagnetic phase (including 
carrier's entropy) is always more stable that the ML po- 
laronic ansatz. Notice that our original statement re- 
mains true: the existence of stable polarons would imply 
thermodynamic instability of the homogeneous phase. It 
only happens that the carrier's entropy makes the para- 
magnetic homogeneous phase to be prefered over the po- 
laronic one at low doping, making the very existence of 
polarons at finite doping uncertain. 

Given the previous situation, we study the existence 
of PS in the pyrochlorcs performing the usual mean field 
calculation without simplifying hypothesis for the carri- 
ers. The Hclmholtz free energy T contains the core Ts 1 



carrier T% 
pressions: 

T s = 



and core-carrier contributions J- S g, with ex- 
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where z is the number of nearest neighbors, T>(e) repre- 
sents the paramagnetic density of states (per spin) and h 
describes the effective field characterizing the mean field 
distribution of core spins, assumed to be classical vectors 
of unit length, therefore, s — ctnh(/i) — h~ l . 

Minimization of T leads to the following equations: 
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with the occupation constraint: 
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We have solved the previous equations for parameters 
expected to apply in the case of pyrochlores: kTo = 4p = 
0.1 t, and J' = t, where t, a measure the electronic energy 
scale, is taken to be the kinetic energy of an electron con- 
fined to a unit cell volume. Searching for negative values 
of /i = d n T, we arrive at the phase diagram plotted in fig. 
H for small carrier concentration. We see that, beyond a 
carrier concentration, a region of negative compressibil- 
ity opens around the ferro-para transition. Therefore, we 
confirm the existence of PS as suggested by the general 
arguments of section II. Notice that in fig. || we show 
the area of intrinsic instability: the standard Maxwell 
construction would increase this PS region even into the 
paramagnetic phase. We have plotted (fig. [j]) the free en- 
ergy of the exact solution (continuous line) , compared to 
the ML polaron ansatz (dashed, straight line) and the so- 
lution ignoring the electron entropy (dashed-dotted line) 
for kT — 0.115 t. As explained at the beginning of this 
section, the correct homogeneous phase (including en- 
tropy) has the lowest free energy even at low concentra- 
tion. In contrast, having ignored the electronic entropy 
would have led to the prediction that polarons are the 
prefered configuration (see two upper curves of fig. ^). 
Nevertheless, PS does exists, and the onset of negative 
curvature (hardly visible to the bare eye) is marked with 
an arrow in that figure. 

For the parameters corresponding to fig. ^, electrons 
are strongly non degenerate, and the numerical results 
almost coincide with the classical limit for the carri- 
ers, whose Ginsburg-Landau analytic treatment we now 
present for completeness. Assuming the electrons to be 
a classical ideal gas and performing an expansion in the 
core s, and electron magnetization s, one obtains the 
following expressions for the contributions to the free en- 
ergy: 
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where ?t,q is a reference density marking the onset of the 
quantum regime. 

Minimization of the free energy leads to the following 
transition temperature: 



kT c = 
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In the vicinity of T c , core and electron magnetization are 
related by 
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and the scaled inverse compressibility at T 
by: 



(38) 



T c is given 
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where k 1 = kT c /n, is the inverse compressibility of the 
paramagnetic phase, k becomes negative beyond nps — 
54T C 4 /(5J' 4 ) = 1.266 x 1CT 3 , for the parameters used 
here, very close to the exact results of fig. |2| 

Aa a final remark, it is important to realize that 
the precise configuration of the coexisting phases in the 
real material will be complicated by effects beyond the 
present treatment (disorder, domain boundary contribu- 
tions ...). Nevertheless, the associated spatial inhomo- 
geneity will certainly affect the physics of these com- 
pounds in an important way. This strongly suggests that, 
even if the original polaron argument is doubtful, the as- 
sociated physics explored by ML might well apply to the 
pyrochlores. 



V. DOUBLE EXCHANGE SYSTEMS. 

The double exchange (DE) model describes systems 
with local spins and itinerant electrons, in the limit where 
the Hund coupling between the electrons and the spins 
is much larger than other scales |34[ . The electrons are 
always polarized in the direction of the local spins, and 
hopping (t) to neighboring ions is modulated by the rel- 
ative orientation of these spins, leading to the following 
Hamiltonian 



n = -t 
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z ia Z ja c l c j 



H.c. 



(40) 



where z ia is the spinor describing the orientation of the 



core at site i: = cos^flj) 



•tj 



sin(|( 



exp(-i^). 

The model, put forward by Zener j54j, has received 
great attention in recent times as the building block 
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FIG. 2. Phase diagrams for the mean field theory of eqns. 
( [jo|]3l| , |32| ) for the parameters explained in the text. The solid 
line is the (scaled) critical temperature versus carrier concen- 
tration. The region of negative compressibility is bound by 
the dashed line and the critical temperature. 



needed to explain the physics of Mn perovskites. We will 
restrict ourselves to its simplest one-band version, leaving 
aside further features such as orbital degeneracy, antifcr- 
romagnetic couplings or Jahn- Teller distortion, probably 
required for a more realistic description of manganites. 

PS in double exchange models has been extensively 
studied It was originally observed in numerical 

studies of the DE model with AF couplings between core 
spins, raising doubts about its existence in the bare ver- 
sion |3||||]]. However, we have recently shown that PS 
tendency is an intrinsic feature of the simplest DE model 
without additional terms 0. For completeness, we in- 
clude here an analysis of the compressibility using the 
Ginzburg-Landau expansion outlined in the preceding 
sections. The approach is a direct extension of that used 
in J35|. We leave out the study of PS at zero tempera- 
ture due to the competition between direct antiferromag- 
netic interactions and the double exchange mechanism 
H , which also can be cast in terms of a Ginzburg-Landau 
expansion. 

The model is determined by the polarization of the 
localized spins, s, the density of states of the itinerant 
electrons in the absence of spin disorder, T>o(e), and the 
temperature T. The free energy of the spins comes from 
the entropy due to their thermal fluctuations only, and 
coincides with eq. ( |34] ) , except for the absence of the di- 
rect exchange term. We assume that the electron gas 
is degenerate, an excellent approximation for tempera- 
tures in the range of T c , and replace its free energy by 
the ground state energy in a background of fluctuating 
spins. The bandwidth of the electrons is reduced by a 
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factor / = (cos(i?ij/2)), where ■dij is the angle between 
neighboring spins. Thus, the electronic energy can be 
written as: 



E = fK = f 



eV Q (e)de 



(41) 



where /iq denotes the chemical potential in the absence 
of spin disorder, and Wo is the lower band edge (note 
that —Wo < Ho < Wo)- We can now expand / in terms 
of the magnetization s, to obtain: 



2 2 , 6 , 



(42) 



which completes the required Ginzburg-Landau expan- 
sion. The Curie temperature of the model is given by 
T c = —jrKo (Kq is negative). The magnetization is: 







\kT+4-K 



rkT- 



which leads to: 
T = 



-k 



Tf(^) T<T C 



T>T C 
T<T B 



(43) 



(44) 



where /(it) = g ^ios^ 

Finally, the scaled inverse compressibility is: 

T -> T c + 
T c -0 



h' 1 



2 

3X> («)) 



14pp rp 



The compressibility below T c is negative when: 

J*eL< V 



(45) 



(46) 



This result is independent of the initial bandwidth, 2Wj. 
In general, there is always PS at low fillings, because 
Mo -> "W while -> 0. At half filling, Mo -» 0, 

and there is no PS. 
If we use: 



2 0^[= 



we find that there is PS if: 



ImoI > Wo 



(47) 



(48) 



which corresponds to a number of electrons per site ~ 0.2. 

The previous expansion in powers of the order param- 
eter is not necessary, and the complete mean-field equa- 
tions can be written as follows 



1 If 1 

x- - = - J djV (j) tanh 



26 



fQ 



c?7 7 T>o (7) tanh 
1 



a + fl 
26 

q + h 

26 



(49) 
(50) 





FIG. 3. Phase diagrams (upper panel) for the mean field 
theory of eqns. (^,^), and (lower panel) for that of eqns. 
( pl|j5^ ), both computed using an elliptic density of states. 
The solid line is the (scaled) critical temperature versus hole 
concentration. The dashed line marks the boundary of the 
region of negative compressibility. 



where s — ctnh (Q) — Q 1 is the magnetization, x 
describes the doping level (electron/hole), and / = 

\J \{]~ + s 2 ) = {cos 2 (dij)) 1 / 2 . The energy, temperature 
and chemical potential are scaled in units of the half- 
bandwidth in the absence of spin disorder by 7 = e/Wb, 
a = fjt/Wo, and 6 = kT/W . 

If T>o is taken to be the density of states of the original 
crystal lattice, neglect of the Berry phase [S7J collected by 
the hopping electron is implicit. More appropriately, we 
can assume that spin disorder cancellations favor retraced 
paths p6|, leading to a Bethe lattice density of states, 
with eq. (|47|) as the infinite coordination limit. In fig. || 
(upper panel) we show the phase diagram (scaled tem- 
perature versus hole concentration x) for the density of 
eq.(^7j). Notice that the PS region corresponds to the in- 
trinsic instability (negative compressibility). Global sta- 
bility (Maxwell construction) would increase this region. 

The simple DE model offers a clear example of the 
general mechanism for PS in the presence of an ordering 
process. In fact, the mean field treatment is intrinsically 
unstable versus PS in the limit of low carrier density. To 
see this, notice that the energy scale of the magnetic or- 
der is the electronic energy, proportional to the carrier 
concentration in the dilute limit. Doubling, for instance, 
the carrier concentration close to the Curie temperature 
produces a finite increase of the bandwidth. The chemi- 
cal potential instead remains tied to the band edge and 
follows its fate: it decreases with increasing electron den- 
sity (thermodynamic instability). In our electron-hole 
symmetric model, the same would apply to the dilute 
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hole limit. 

Although we do not expect the mean field approxi- 
mation to be a serious drawback, certain limitations are 
clear. For instance, the inability to distinguish between 
short-range correlations and long-range order confines 
the negative compressibility region to the ferromagnetic 
phase, for only there does the bandwidth change with 
magnetic order. In what follows, we will show that PS 
is a robust feature that survives several improvements of 
this basic mean field approach. 



We have solved the mean field equations for a semi- 
elliptic density of states and S — 3/2 M. Several fea- 
tures of the solution are shown in fig. |J| (lower panel). 
While the critical temperature is hardly affected, the re- 
gion of negative compressibility enters now into the para- 
magnetic phase, an expected physical feature not present 
in the previous approximation. In this approach, the 
fermion bandwidth collapses to zero at high temperature. 
Although this is likely an artifact of the approximations, 
it takes place well above the Curie temperatures to be a 
source of concern. 



A. Schwinger bosons. 

An improved mean field theory can be implemented 
using the Schwinger boson |38|,f7| method. A nondy- 
namical field Ai enforces the constraint z ia z i(T = 1 at 
every site. The core spins are quantized according to 
[z ia .,Zj a ,] — SijSa-a-' /2S. The mean field Hamiltonian is 
obtained through a Hartree decoupling of the bosonic 



z i<j z jo an( i fermionic qCj hopping terms: 



(51) 



(ij) {ij),<r 



where K — (cjc.) and / = (z ia Zj a ) Such a model was in- 
troduced by Sarker j39j , who identified a Curie transition 
and found that the e g fermion band becomes incoherent 
above T c . Accounting for the possibility of condensation 
of Schwinger bosons, we write = { z u )■ Assuming 



condensation only at k — 0, we define p 
The mean field equations are then 



l**=o l0 -l 



= P 



f = P + 



^/W)ctnh(^ 

If 1 . . /A - if 7 

— / c?7 7X>(7)ctnh 



4S*6 



1 

" 2 

K 



(ca -f- fry 

i J d 11 V( 1 ) tanh/ 



v 29 

a + fl 
26 



(52) 



where 7 = e/W , A = A/W , 6 = k B T/W , and a = 
fx/Wo, Wq being the half-bandwidth in the absence of 
spin disorder. 

Two aspects distinguish this treatment from that of the 
previous section. First, core spins are quantum objects 
with intrinsic dynamics: a minute correction close to the 
Curie temperature. Second and more important, long- 
range magnetic order (condensation of Schwinger bosons) 
is not a requisite for short range correlations: the latter 
exist and change with temperature even in the paramag- 
netic phase. 



B. Critical spin fluctuations 

The rapid change of spin correlations ties the PS re- 
gion to the magnetic ordering transition. One can worry 
about the importance of critical spin fluctuations, ne- 
glected in the molecular field approximation. Assuming 
a nearest- neighbors, tight-binding lattice with hopping 



amplitude U 



-t cos(0y/2), and expanding the half- 



angle between spins in the usual manner: cos(#/2) ~ 
a + axcos(0), the standard virtual crystal decoupling of 
the electron-spin system allows integration of charge de- 
grees of freedom . This produces an effective Heisenberg 
model for the core spins (assumed classical for simplicity) 
and the total free energy (charge and core spins) is given 
by 



F = K^x) + F H {J,T) 



(53) 



where Koo is the electron energy (degenerate limit as- 
sumed) of the fully disordered limit tij = a^t. !Fh de- 
scribes the free energy of the classical Heisenberg model, 
with coupling constant J = 2ai|(cjc,,}|i, where tilded 
quantities are electron properties scaled to unit half- 
bandwidth. A lengthy but simple calculation permits 
de following connection between the stability criterium 
and thermal properties of the Heisenberg model (cubic 
lattice): 



1 



3 2 T 



dn 2 T^{p) 



with 



and 



9{x) 



[ 1 - g(x) f H (T) ; 



2aijl 2 V(il) 



f H {T) = 



6|(4 Cj )| 



TC(T) 



J [a Q + ai(cos(%)) 



(54) 



(55) 



(56) 



fi is the chemical potential and x the carrier density. 
C(T) represents the specific heat of the Heisenberg 
model. The expected increase of C(T) near the criti- 
cal temperature T c links the tendency towards instability 
with magnetic ordering. In any case, for a fixed value of 
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T/T c , the system is bound to exhibit PS, (dfx/dx) < 0, 
owing to the diverging behavior of g(x — > 0). Extract- 
ing the temperature dependence of fji (T) from published 
Monte Carlo data (40) for the cubic lattice, we estimate 
PS accompanying the ordering transition up to x ~ 0.11. 
The resulting phase diagram and intrinsic instability re- 
gion are shown in fig. || (upper panel). Notice that no ap- 
proximate treatment is assumed for the criticality of the 
effective magnetic Hamiltonian, the only (major) surviv- 
ing approximation being the virtual crystal (mean field) 
decoupling of mobile carriers and core spins. 



C. Infinite dimension 

The mean field decoupling of charge and core spins 
fluctuations remains in all the previous analysis. This 
manifests itself in a mere rescaling of the electronic den- 
sity of states with magnetic order. As the PS instability 
can be seen as a manifestation of correlation between 
charge and spin at a macroscopic level (the only possi- 
bility left in a mean field approach) concern might arise 
about its fate when charge-spin coupling is properly con- 
sidered. Although in the general case a numerical ap- 
proach would be necessary to answer this question, there 
is a limit in which the problem is solvable but the charge- 
spin coupling remains non trivial: infinite coordination 
(dimension). The relevance of this limit for the under- 
standing of correlated systems in 3-D has been recognized 
in recent years [ p]j42| . 

In the DE case, this problem has been considered by 
Furukawa |4^ |, and indeed PS has been observed In 
those works, the DE model included a large though finite 
Hund coupling, leading to AF correlations between core 
spins that certainly increase the tendency towards PS. 
This fact, together with accummulated experience from 
numerical studies of doped antiferromagnets, has led to 
the extended belief that antiferromagnetism is required 
for PS to exist. We now show that, in spite of recent 
suggestions on the contrary, the simple DE model with- 
out AF additions indeed exhibits PS close to the ordering 
temperature, according to the general explained in this 
paper. 

The average density of states for electrons moving in 
an infinitely coordinated Bethe lattice according to the 
DE Hamiltonian is given by: 



(X>(e)) fi = -(l/7r)3<fl( e> n))tf 
where the local Green's function satisfies: 

g{e,(i)= [e-£M) 
and the self-energy is determined by: 
1 + • ft' 



(57) 



(58) 



£(e,fi) = (Wo/2) 2 ( 



e-E(e,fi') ) 6 





0.00 0.02 



FIG. 4. Phase diagrams (upper panel) for the cubic lattice 
in the virtual crystal approach, including MonteCarlo data 
for the Heisenbeg model (subsection B) and (lower panel) for 
the Bethe lattice of infinite coordination. The solid line is 
the (scaled) critical temperature versus hole concentration. 
The dashed line marks the boundary of the region of negative 
compressibility. 



is the unit vector describing the local core (classi- 
cal) spin, Wq represents the half-bandwidth for the fully 
aligned case, and angular averages are taken with the 
probability distribution P(f2). Notice that, though the 
probability distribution for spins in different sites factor- 
izes (confining PS the the ferromagnetic phase), charge 
and spin fluctuations remain coupled at the same site, 
and the density of states is not merely rescaled by mag- 
netic order (eq. |59| ). 

The complete thermodynamic problem is solved by 
finding the probability distribution that minimizes the 
total free energy. 



T = T e] - TS 



(60) 



(59) 



where T e i is the free energy of fermions for the average 
density of states and S is the core spin entropy, both 
depending on V(Vl). 

For computational convenience, we have solved the 
problem with two additional features: degenerate limit 
for carriers and a core spin distribution parametrized 
by an effective magnetic field (detailed analysis of the 
exact solution for a few temperatures shows the previ- 
ous simplifications to introduce only minute corrections) . 
The generic results of previous treatments are reproduced 
here. This is shown in fig. || (lower panel), where the 
phase diagram presents the expected region of intrinsic 
instability (negative compressibility) at low carrier den- 
sity. The inclusion of a finite Hund coupling will cer- 
tainly increase the PS region, as observed in previous 
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kT/ Wo = 0.01 5 



it 





1 1 1 1 1 1 ' 1 

0.0 0.1 0.2 0.3 

X 

FIG. 5. Compressibility of the DE model for the (infinitely 
coordinated) Bethe lattice, in units of the compressibility of 
a free fermion system with identical density of states, versus 
hole concentration at constant temperature. 



o ' ' — ' ' — ^ — ' ' ' — ' ' ' — ^ — ' — ' ' — ' ' ' — ' ' ' 

O 0.005 0.01 0.015 0.02 

T / W 

FIG. 6. Region of phase separation, as function of mag- 
netic field, in the double exchange model, for an electron den- 
sity, x — 0.10. Temperature and magnetic field are normal- 
ized to the electron bandwidth. 



works |£]]6|]. All these studies support the view that the 
general scheme for PS presented in this paper applies to 
the simple DE model with features which are robust and 
not an artifact of approximations. 

As mentioned before, PS can be thought of as a long 
range instability caused by the coupling between carriers 
and core spins. Such mechanism is operative even in the 
region where the system does not phase separate, and the 
carriers compressibility can be much enhanced close to 
the Curie temperature. This statement does not contra- 
dict the fact that the density of states at the Fermi level is 
featureless and shows no major change around the Curie 
temperature jf4|. Remember that, if the bandwidth were 
not affected by the spin order, the compressibility would 
be merely proportional to the density of states at the 
Fermi level, that means featureless. It is the bandwidth 
change with magnetic order what produces this enhance- 
ment in the region close to the onset of magnetic order. 
This is shown in Fig where the compressibility of the 
infinitely coordinated Bethe lattice is measured in units 
of the compressibility of a free fermion system (that is, 
no carrier-spin coupling) with the same density of states. 
This behavior is not a peculiarity of the infinitely di- 
mension limit, and remains the the same in all previous 
approaches. The possible connection between this en- 
hancement and the unusual properties of the manganites 
remains an open and intriguing question [^) . 



temperature dependence. 

The increase in the compressibility which leads to 
phase separation is asociated with the strong coupling be- 
tween magnetic and charge fluctuations, as measured by 
the dependence of the critical temperature on electronic 
concentration. The magnetic field suppresses magnetic 
fluctuations leading to a smaller increase in the compress- 
ibility near T c . Phase separation should dissappear at 
sufficiently large fields, when the temperature no longer 
induces significant magnetic fluctuations. 

A generic case which shows the dependence of the re- 
gion where phase separation occurs on magnetic field, at 
fixed nominal electronic concentration, is shown in fig. 
(g). The calculations have been done for the double ex- 
change model of the previous sections, using the mean 
field equations ( P4 fl3| ), plus an applied field. The region 
of negative compressibility is shown. As in previous ex- 
amples, a Maxwell construction gives a somewhat larger 
region. 

In a small applied field, phase separation takes place 
above the Curie temperature. This is due to the round- 
ing of the discontinuity at T c of the compressibility in- 
duced by the field. The region of phase separation lies 
between a high magnetization and a low magnetization 
phase, in close analogy with the phase diagram of an or- 
dinary liquid-vapor transition. At high fields, phase sep- 
aration is completely suppressed. At the highest possible 
field we find a critical point. 



VI. INFLUENCE OF A MAGNETIC FIELD 

So far, we have considered the existence of phase sepa- 
ration near a ferromagnetic (or antiferromagnetic) - para- 
magnetic transition, at zero magnetic field. When a field 
is applied, the features associated to the transition are 
smoothed, and, at sufficiently large fields, the magnetic 
moments are aligned, and the magnetization has a weak 



VII. PHASE SEPARATION AND DOMAIN 
FORMATION 

The analysis in the previous sections suggests that PS 
should be a frequent feature of magnetic transitions. PS 
on a macroscopic scale, however, cannot occur, as it re- 
quires an infinite amount of electrostatic energy. We can 
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extend the mean field analysis used in the previous sec- 
tion to study the role of the Coulomb interaction. If we 
first neglect it, the charge susceptibility obeys: 



lim Xo(q) 

q^O 



-V 



,\dn 

9yU 



(61) 



where V is the volume of the unit cell. A charge 
fluctuation p(q) induces an electrostatic potential, 
(4-7re 2 /q 2 )p{q), which induces more charge polarization. 
The RPA selfconsistent equation for the charge polariz- 
ability becomes: 



lim Xrq) = — 
g^o {HJ 1 



Xo(«) 



-v 



-1 2 

n k 



47re 2 

q' 2 



X (q) 



1 



47re 2 
4 1 



v- 



(62) 



when, in the absence of Coulomb repulsion the system 
exhibits the instability, k < 0, the denominator has a 
pole at: 



47re 2 



(63) 



and charge fluctuations of shorter wavelengths remain 
unstable, whereas the long range nature of the Coulomb 
term prevents the formation of larger charge inhomo- 
geneities. Thus, we expect the formation of domains at 
length scales comparable to q~ x . As mentioned previ- 
ously, this analysis does not take into account neither the 
cost in magnetic energy associated with the formation of 
domain walls nor the effect of impurities. In addition 
to macroscopic charge neutrality, these are likely to be 
among the major factors affecting the spatial coexistence 
pattern in real materials. 



VIII. CONCLUSIONS 

We have discussed a general framework which shows 
that PS is likely to occur near magnetic phase transi- 
tions. The dependence of the critical temperature, or 
critical couplings, on electronic density leads to a reduc- 
tion to the compressibility in the ordered phase. This re- 
duction tends to be of the same order as the value of the 
compressibility in the disordered phase. The existence 
of PS typically depends on numerical constants of order 
unity, related to the electronic structure. The coupling 
constants do not require a special fine tuning for PS to 
take place. The analysis presented here is probably rele- 
vant to understand a variety of experimental findings in 
the manganites, the pyrochlores and doped antifcromag- 
nets. 

We have studied PS mostly within a conventional mean 
field framework. Corrections due to critical fluctuations 
will usually enhance the effects discussed here. The com- 
pressibility is associated with the second derivative of the 
free energy with respect to electronic density, in a simi- 
lar way to the specific heat. Critical fluctuations tend to 



increase the divergences of the derivatives of the free en- 
ergy at the transition, so that the estimates reported here 
probably underestimate the tendency towards PS. The 
analysis also shows that, even in the absence of phase 
separation, an enhanced coupling between charge and 
magnetic fluctuations is expected, due to the reduction 
in the charge compressibility. 

Finally, we have discussed the way in which electro- 
static effects frustrate the formation of macroscopic do- 
mains, and we give a scheme to calculate the scales at 
which domain formation is expected. 
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